Existence of vertical spin stiffness in Landau-Lifshitz-Gilbert equation in 

ferromagnetic semiconductors 
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We calculate the magnetization torque due to the spin polarization of the itinerant electrons by 
deriving the kinetic spin Bloch equations based on the s-d model. We find that the first-order 
gradient of the magnetization inhomogeneity gives rise to the current-induced torques, which are 
consistent to the previous works. At the second-order gradient, we find an effective magnetic field 
perpendicular to the spin stiffness filed. This field is proportional to the nonadiabatic parameter /3. 
We show that this vertical spin stiffness term can significantly modify the domain-wall structure in 
ferromagnetic semiconductors and hence should be included in the Landau-Lifshitz-Gilbert equation 
in studying the magnetization dynamics. 

PACS numbers: 75.60.Ch, 72.25.Dc, 75.30.Gw, 75.50.Pp 
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I. INTRODUCTION 

Ferromagnetic systems have attracted much interest 
for a long history because of the intriguing physics and 
applications, 1,2 As the development of information tech- 
nology, the research on magnetization dynamics in mi- 
cromagnets has become an active field Great efforts 
have been devoted to this field by aiming to manipulate 
magnetization more efficiently^— For theoretical simula- 
tion, the magnetization dynamics is usually described by 
the Landau-Lifshitz-Gilbert (LLG) equation;^ 

ri = —711 x H c ff + an x h — (1 — /3nx)(v s • V)n, (1) 

where n represents the direction of the magnetization. 
H e g in the first term on the right-hand side of Eq. (JTJ 
is the effective magnetic field which drives the magneti- 
zation procession and determines the domain structure 
in the equilibrium states. Up to date, different sources 
of effective magnetic field have been identified, e.g., the 
external magnetic field, the crystal anisotropy induced 
by the spin-orbit coupling, the exchange energy due to 
the spatial inhomogeneity, and also the demagnetization 
field from dipole-dipole interaction^ The second term, 
the Gilbert damping torque, describes the magnetization 
relaxation to the effective field axis on the time scale of 
l/(ajH c s)Jt The last one with first-order gradient of the 
magnetization is the spin torque induced by the trans- 
port of the spin polarized itinerant electrons, where v s 
is proportional to the spin current density*^— The fi- 
terni ) 12 ' 13 first proposed by Zhang and Lijis was demon- 
strated to be critically important to the current-driven 
domain wall motion, which overwhelms the threshold 
current due to pinning force and transverse anisotropy 
for domain wall motio n 10 i 18 ' 19 and results in the steady 
domain wall velocity cx /3/a in the absence of the exter- 
nal magnetic field. Therefore, the determination of the 
effective magnetic field and the parameters in the LLG 
equation, such as a and /3, is an important issue for mag- 
netization dynamics studyi 12 ' 14 ' 15 ' 20 ' 21 



Previously, we have derived the Gilbert damping (a) 
term based on the kinetic spin Bloch equation (KSBE) 
approach 2 ^ in homogeneous ferromagnetic systems ^ In 
the present work, we derive the whole LLG equation from 
the s-d model in inhomogeneous ferromagnetic systems 
based on the same approach. From the first-order gra- 
dient of the magnetization inhomogeneity, we obtain the 
current-induced torques which are consistent to the above 
LLG equation. Within the second-order gradient, we find 
that the LLG equation should be modified and written 
as 

ri = —711 x H c ff + an x ri — (1 — /3nx)(v s • V)n 

-^A ss nx(l-/3nx)V 2 n. (2) 

The second-order gradient introduces two contributions 
to the effective magnetic field. The one given by A ss V 2 n 
is identified as the spin stiffness field discussed in pre- 
vious works . 10 ' 25 ' 26 The other term, in the form of 
— f3A ss n x V 2 n, has never be referred in the literature. 
In the present paper, we call it "vertical spin stiffness" in 
the sense of the fact that this new field is vertical to the 
plane defined by the magnetization and the normal spin 
stiffness V 2 n. Interestingly, this vertical spin stiffness 
can not be written in terms of the free energy, and there- 
fore, it can not be derived from the functional derivative 
of the free energy with respect to the local magnetiza- 
tion H cff = -6F[M d ]/M d previously^ We find that 
this vertical spin stiffness results in the tilt of the mag- 
netization. The new term can significantly change the 
domain-wall structure in ferromagnetic semiconductors. 
Since the magnitude of this field is proportional to the 
factor of /3, the proposed effect is expected to be impor- 
tant in ferromagnetic semiconductors where /3 is large 
due to the strong spin-orbit interaction! 20 ' 23 

This paper is organized as follows: In Sec. II, we setup 
our model and derive the KSBEs for the itinerant elec- 
trons in the inhomogeneous ferromagnetic system. We 
calculate the spin torque by solving the KSBEs in Sec. Ill 
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and discuss the results in Sec. IV. Finally, we briefly sum- 
marize in Sec. V. 



II. KSBES 

We use the exchange interaction Hamiltonian density 
-Hsd(r) = Mn • cr with M denoting the coupling con- 
stant. Following Ref.Hil we here assume that the ferro- 
magnetic interaction exists among n(r) and show later 
that this assumption is justified after integrating out the 
conduction electrons (giving rise to spin stiffness). The 
Pauli matrices cr are used to describe the itinerant elec- 
trons. In contrast to the previous work on the Gilbert 
damping^ we introduce the inhomogeneity by consid- 
ering the position dependence of the magnetization di- 
rection n(r) = M,j(r)/M|i. Md is the uniform satu- 
rate magnetization. For the strong exchange coupling 
in ferromagnetic system, the rotation framework^^ is 
employed here. That is, the local spinor operators of the 
itinerant electrons are defined as a(r) = (a-)-(r), aj_(r)) T , 
with f (i) labeling the spin orientation parallel (anti- 
parallel) to n(r). Therefore, one has H s & — Ma*a z a. 
The spinor operators a(r) are connected to the ones 
defined in the lattice coordinate system c = (c^,cjJ T 
via the unitary transformation a(r) = U(r)c. The 
transformation matrices is given by U (r) = m(r) • cr 
with m = (sin(0/2)cos<^, sin(0/2)sinyj, cos0) for n(r) = 
(sin^cos^, sin^sin^, cos9)^^ 

In the rotation framework, the kinetic Hamiltonian 
should be expressed as = |Vc| 2 /2m = |(V + 

iA)a\ 2 /2m = |Va| 2 /2m + Ha, where the gauge field in- 
troduced by the coordinate transformation is given by 
Ai — —iWdiU — (m x Sfm)j<7j = A\oi. Then, one ob- 
tains the Hamiltonian density associated with the gauge 
fiekU 



i 

2m 



^[aU;Via- (Via)Uia] 



A 2 



a ' - — a. (3) 
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To derive the KSBEs of the itinerant electrons, we fol- 
low the nonequilibrium Green function approach . 22 ' 27 ! 28 
The Dyson equation of the contour-ordered Green func- 
tion can be expressed as 



G(l,2) = G (l,2) + J d3G (l,3)t/ e (3)G(3,2) 
+ J d3 Jd4G (l, 3)S(3, 4)G(4, 2) 
= G (l,2)+^d3G(l,3)[/(3)G (3,2) 

d3 /cMG(l,3)E(3,4)Go(4,2), (4) 



c •/ c 



where the single particle contour-ordered Green func- 
tion between two space-time points (1) = (ri,ii) and 
(2) = (r2,*2) on the contour G is defined as G(l,2) = 



— i(Tc['0_ff(l)V-'^(2)])- 27 U e describes the local electric po- 
tential energy, whereas E represents the self-energy cor- 
rection due to interactions, such as the electron-impurity, 
electron-phonon, and electron-electron interactions. Go 
stands for the free-particle Green function. The left- and 
right-inverses of Go(l, 2) are given by 

(G (xx,x 2 ))zl = id tl -H Q {vx,rx), (5) 

(Go {xx,x 2 ))-, = -i dt 2 -H (- P2,r 2 ), (6) 

with Hq = + H s< x- By multiplying them to Eq. ((4]), 
one obtains 

i(d fl +d t2 )G(l,2) = [ff (Pi,ri) + [/ e (l)]G(l,2) 
-G(l,2)[# (- P 2 ,r 2 ) + C/ e (2)] 

d3[£(l,3)G(3,2)-G(l,3)£(3,2)], (7) 

with pi corresponding to the momentum operators —id ri . 
To transform the above equation to the center-of-mass 
and relative variables 



R = (ri + r 2 )/2; r = n r 2 ; 
T=(tx+t 2 )/2; t = tx-t 2 , 



(8) 



one rewrites the following Green functions as 

G(l,2) = G(R,r,T+f,T-f), (9) 

G(l,3) = eT 9 -G(R,r 1 -r 3 ,r+|,i 3 ), (10) 

G(3,2) = e^ 9 -G(R,r3-r 2 ,i3,r-|). (11) 

The self-energy can be written in the same way. Simi- 
larly, one obtains 



U e (l) = e^ R £/ e (R,T+f), 
U e (2) = e-^[/ e (R,T-§). 

The Hamiltonian can be written as 

flo(Pi,ri) = e^°H o aP R + p,K), 



(12) 
(13) 

(14) 



-ffo(-P 2 ,r 2 ) = ^o(4PR + P. R ) e " § R - (15) 

where Pr and p represent the momentum operators re- 
spect to the center-of-mass and relative variables. Ac- 
cording to Eq. , the left- and right-operators of the 
kinetic Hamiltonian should be different and can be writ- 
ten as 

ff k (p,r) = ^[p 2 -*(V-A) + 2A-p + A 2 ],(16) 

if k (-P,r) = ^[p 2 +i (V-A)-2p -A + A 2 ](17) 

separately. By assuming the magnetization varies 
smoothly respect to the spatial coordinates, we only keep 
the spatial gradient up to the second order. Since the 
gauge field is already the first-order gradient, both V ■ A 
and A 2 are second-order gradient terms. We include both 
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of them and neglect the higher-order ones. The gradient 
term of the gauge field can be written in the center-of- 
mass coordinate system as 

V ri -A(n) = ±V R - A(R+§) + V r -A(R+§) 
= iV R -A(R) + d ri ^.^(R) 
= Vr-A(R). (18) 

Similar calculation gives 

V r2 • A(r 2 ) = Vr-A(R). (19) 

Moreover, one can easily show that A 2 (r;) w A 2 (R) — 
IJ^i i(A\) 2 with I representing the unit matrix. 

By substituting all these equations into Eq. ([7|) and 
doing Fourier transformation respect to the relative co- 
ordinate r, one obtains 

tftrG(R,k,ii > t a ) 

= e l ^ dkd ^"~ d ^ ( ">Ho(k, R)G(R, k, t u t 2 ) 

- e -<i(a.a£° - 3 R 9 f°)G(R, k, * x , t 2 )# (k, R) 
+ e 4 W Ue {K, ti)G(R, k, h,t 2 ) 

_ e -4^^ e G(R, k, tx, h)U e (K, t 2 ) 

+ J *3[e^ ( ^ 9l ~ 9fag) S(R, k, h,t 3 )G(R, k, t 3 ,t 2 ) 

- e^ a £-«>G(R, k, h, * 3 )E(R, k, t 3 , t 2 )]. (20) 

The details can be found in Appendix A. We then per- 
form the gradient expansion up to the first order and 
obtain 

id T G = [H + U e ,G}~ H^Ho, d R G} + f{d R d k ff , G} 
+ |{5r(^ + U e ), d k G} + Jdt 3 (XG - G£){21) 

where all the quantities are defined at R and k. We 
should point out that the commutator notation [Hq + 
U e , G] is still used although the left- and right-operators 
of Hq are in different expressions [see Eqs. (16) and (17)]. 
By taking the isochronous condition, i.e., t — > 0, one 
has22 

J dt 3 [Z(T, r)G(r, T) - G(T, r)E(r, T)]< 

= f dr[E>(r,T)G < (r,r)-S<(r,T)G > (r,r) 

- G> (T, r)E< (r, T) + G< (T, r)E> (r, T)] , (22) 

where the lesser and greater Green functions are de- 
fined by G<(t,t') = *(^(t / )^(t)) and G>(r,r') = 
—i(ipH(T)ijj H (T')). Therefore, the correlation function 
G < (T, T) can be written as 

id T G< = [ dr{T, > G < -T, < G > -G > Y, < +G < Yr > ) 

J —OO 

+ I {Or^Ho, G<} + 1 {a R (i/ + c/ e ), 9 k G<} 

+ [H + U e , G<] - ^{dkHo, 9rG < }. (23) 



Within the generalized Kadanoff-Baym ansatz^Z. we have 
G<(R,k,T,T) = ip k (R,T) where j0 k (R,T) is the local 
density matrix of the itinerant electrons with momentum 
k located at R. Therefore, we write the general form of 
the KSBEs of the itinerant electrons as 

d t Pk + i[H ,Pk] + H v k#0; V Rj o k } - V R U e ■ V k p k 

- |{V R ff , V ki o k } - ±{V R • V k iJ , Pk} 

= d t pk[ c scat + 9tPk|Lt- (24) 

On the right-hand side of above equations, <9t/9 k |g Cat and 
c^p k |s cat from the integral term in Eq. ([23| represent the 
spin-conserving and spin-flip scatterings. The details of 
these terms can be found in Ref.l22l 

We specify the Hamiltonian H = Mu z + [(k 2 + 2A t ki + 
A\A\)I^f i(V ■ A)] /(2m) with the upper (lower) sign rep- 
resenting the left (right) operator case. The electric po- 
tential energy is given by U e = eE • R (e > 0). The final 
form of the KSBEs is given by 

d t pu + i[M<r g , p k ] + [Ai,pk] + 2^{A, V R( o k } 
+ ^k ■ V R p k - eE • V k( o k - !HV k p k , VrA,} 
= ftp k | s c cat +%>k|Lt- ( 25 ) 

Interestingly, we find that the contribution from — ^{V r • 
V k iJo, Pk} is completely canceled by the gauge field gra- 
dient term from i[H , p k ] and the A 2 -term is irrelevant. 

III. SOLUTION OF KSBES 

In general cases, the KSBEs are too complicated to 
solve analytically and the numerical scheme should be 
employed. However, the analytical solution can still be 
expected within some simplification of the KSBEs. In the 
following, we assume: (i) the spatial dependence is weak 
in the rotation coordinate systems, hence VRp k ss 0; 
(ii) the scattering is strong enough to set up the steady- 
state condition <9tp k « 0. Without loss of generality, we 
take the magnetization gradient along arbitrary direc- 
tion. The external electric field is applied on the purpose 
of producing current-induced magnetization dynamics. 
Therefore, the KSBEs become 

-eE ■ V k p k + i[^A t + Ma z ,pk] - |^{V kj o k , VrAJ 

= %>k|scat + dtPklLxt- (26) 

For the steady-state situation with a small current due 
to a static electric field, we assume that the distribution 
of the itinerant electrons is not far away from the Fermi 
distribution. The scattering effect of the spin-conserving 
process is introduced by the relaxation time approxi- 
mation with the average momentum relaxation time r. 
Therefore, one can linearly expand the density matri- 
ces by considering the drift effect, p k — p^I + S k • er = 
Pk + e7 ~E ' ^kP k , where the isotropic density matrices 
Pk = Pk°^ + ^fe ' °" representing the spin polarized Fermi 



4 



distribution in the absence of the electric field. By sub- 
stituting these density matrices into Eq. (j2l))) . the driving 
term and the spin-conserving scattering term cancel out. 
Further, one introduces the average spin relaxation time 
t s and rewrite the spin-flip scattering as 



(s k 



Here, S k describes the equilibrium spin polarization due 
to the spin-splitted band structure in the ferromagnetic 
system. 

Under the above procedures, one finally obtains the 
equations of the steady-state spin polarization 



(2M+^A i )xS k -^V k/ ^-V R A i = - 



> (27) 



Since the s-d exchange interaction can be equivalently 
written as H s d = M-(er) = ^M d -S, the spin procession 
field induced by the s-d exchange interaction — 2MSi/Md 



H X = 5 ^(l + /? 2 )- 1 [^(l-/?nx)V 2 n 

-£(/3 + nx)(j.V)n], (33) 

with (3 = l/(2Mr,). 

We finally obtain the LLG equation with Gilbert 
damping torque as Eq. ((2]), where v s = jP/[2eSd(l + /3 2 )] 
and 



by using the relation { V k /9 k , Vr^L} = 2V k /O k • 
VR^fa/j + 2V k S£ ■ V R Af. Here, we denotes A, = 

Then, the equation of the total spin polarization S can 
be obtained by summing Eq. (|27[) over the k-space 



A ss = n/[4m(l + /? 2 )]. 



(34) 



Here, the effective magnetic field H e g includes the 
isotropic and demagnetization sources 2 ^ as well as the 
external magnetic field. 



-2MxS-|A, x (2^ k fciS k ) 



— s-s 



^F + ^E^V.^-VrA, 



(28) 



For the lowest order approximation, one substitutes S k = 
Sj.+erE- V k S k into Eq. {25]). By considering £ k kS l k = 
0, one obtains 

S-2r s MxS+2r s Aix(ier^S) = S^r^^A,, (29) 

in which the relation ^ k kid^gCk.) = — Sij ^ k g(k) is 
used. The quantity n = 2 ^ k p k stands for the density 
of the itinerant electrons. The spin polarization can then 
be found 



S = 



y + Axy+(A-y)A 
1 + \A\ 2 



(30) 



\Pni with the 



where A = 2r s M - 2r s r^i?jAj and y 
The equilibrium spin polarization S e 
value of the spin polarizability along the magnetization 
direction P is negative, because S e is anti-parallel to 
and M. With the current defined as j = e 2 nE/m, the 
transverse spin polarization is given by 



(l+4r 2 Af 2 )- 1 (-^A^- 
+ x A^ - 



x a \± 

m 1 1 

(31) 



Here, we have neglected the gauge field in the denomi- 
nator. We should point out that the subscript i refers 
to the lattice coordinate axis. In contrast, the gauge 
field is defined in the rotation frame, hence one needs to 
transform the relevant terms back to the lattice coordi- 
nate system. The rotation transformations are given by 
RA^ = —in x d^n and R(z x A ) = \d^n^- One then 
obtains the transverse spin polarization in the lattice co- 
ordinate 



(l + 4r 2 M 2 )- 1 [f^nx V 2 n 



^(j-V)n- 



4m 
MtIp 



Mmt ^2 

° V n 



n x (j • V)n] . 



(32) 



IV. DISCUSSION 

We discuss our results based on Eq. (J2J. One notices 
that the third term on the right-hand side of the equa- 
tion is the current-induced torque obtained in the pre- 
vious workSi 12 i 14 ' 15 The fourth term is independent of 
the current but associated to the second-order gradient 
of the magnetization. This term contains two contribu- 
tions. The one in the form -jgjA ss V 2 n is identified as the 
effective spin stiffness.— In the limit f3 <C 1 or Mr s 3> 1, 
from Eq. (j34[) one has A ss = n/(4m), which is consis- 
tent with the previous result^ However, A ss should be 
modified for finite (3. This stiffness is widely used in 
the study on the domain wal l 10 ' 12 and demonstrated to 
be critical to determine the width of the domain walli 2 ^ 
One finds that the spin stiffness increases with increas- 
ing the exchange coupling strength M, which agrees with 
the previous computation^ By taking n ~ 10 20 cm~ 3 , 
/3 ~ 1 (Ref.[20|), and m — 0.5m e with m e representing 
the free electron mass, one estimates the spin stiffness 
yl ss ~ 1 pj/m in GaMnAs^ However, the other effec- 
tive field in the form —-M-n x V 2 n has not been studied 

Mi 

yet. It is obvious that this torque prevents the magne- 
tization from varying in a plane and induces transverse 
component instead. One notices that this term has no 
contribution to the free energy, because it is always per- 
pendicular to the magnetization. Therefore, it can not 
be derived from the variation of the free energy with re- 
spect to the magnetization, which can explain the reason 
for missing this term in previous works. In the follow- 
ing, we focus on the steady domain-wall solution of the 
LLG equation in the absence of the current to illustrate 
the effect of this new effective magnetic field due to the 
vertical spin stiffness. 

In the ferromagnetic thin film or nanowire structures, 
one takes the total effective magnetic field, H^* = 
Kn x -k - K^n z z + A cff V 2 n - f3'A cS n x V 2 n, l ° with 
K and K± representing the anisotropy constant and 



demagnetization field, respectively. Here, we have 
added stiffness constant ^4g S arising from the non- 
itinerant-electron origin (such as dipole-dipole inter- 
action between the localized d electrons) to describe 
general systems, A e g — (A ss + A® s )/Md and (3' — 
/3A ss /(A e ffMd). The magnetization direction is given by 
n = (cos^sinflcos^jSin^sint/?). Obviously, the ground 
state is the homogeneous configuration with the magneti- 
zation pointing to the easy axis, i.e., the ±x-axis. How- 
ever, the inhomogeneous configurations can also stably 
exist, for example the domain wall structure. To discuss 
the formation of the inhomogeneous magnetization struc- 
ture, we first write the equation of motion for 9 and tp in 
the absence of the current, 

9 + as\\\9tp — — j K ±sm9smtpcostp + d x (sin 2 9d x tp) 
- jP'A eS [d*6 - sm9cos9(d x p>) 2 }, (35) 
— a9 + sm9tp = ryK±sm9 cos 9 sin 2 tp + jK sin 9 cos 9 
- 1 A cS [d 2 x 9-sm9cos9(d x p) 2 } 
- 2 ^ L d x (sm 2 9d x tp). (36) 

For the steady state, one obtains 

d 2 9 = \-jr— sin#cos# + -^sin^sint^cos^sini/? — ft'costp) 
+ (1 + /3' 2 )sim9cos0(c^) 2 ] /{I + f3' 2 ), (37) 

d 2 tp = [^-^-cos# + s'mtp (f3' cos9smtp + cosy) 

- (1 + /3' 2 )2cot9d x 9d x tp] /(l + f3' 2 ). (38) 

At /3' = 0, one obtains a single wall solution in the x-y 



plane (located at x — Xq), tp 



ln(tanf) = 



w 



with W = y/A cS /K. 

When p' ^ 0, the magnetization can not vary in a 
fixed plane since 9 and if are coupled. Unfortunately, 
Eqs. (36) and (37) cannot be solved analytically in gen- 
eral. However, in the absence of the demagnetization 
field (K± = 0), there is a solution where the gradient of 
tp is a constant, i.e., d x tp — A. In this case, the equations 
can be written as 





di9 = 







FIG. 1: (Color Online) The variation of the direction of 
magnetization along the x-axis is plotted as: (a) /3' = 0, 
K± = (or 3K); (b) 0' = 1, K± = 0; and (c) /3' = 1 
and K± = 3K. By denoting the magnetization direction 
n = (cos 6, sin 6 cos tp, sin 6 sin tp) , the position dependences of 



[^sin6>COs6> + A 2 (l + ft )sin0COS0j39) 6 and tp are given in (d)-(f). W = y /A cS /K. 



TT ^ ) [^cos9~X(l+/3' 2 )2cot9d x 9]. (40) 



Obviously, both equations give the solution in the same 
form 



The domain wall thickness is enhanced by j3' as 



ln(tan§) = ^, (41) 

which is just the domain wall solution with the corre- 
sponding width W£ = ^(1 + P' 2 )/[^ + (1 + /3' 2 )A 2 ] 

2|A|^ c ff(l + P' 2 )/{P'K), respectively. The 
Wi determines the value 



w _ l Aclt 2(l+/3' 2 ) 



(43) 



and W£ 

self-consistent condition W% 
of A as 



(42) 



This solution indicates that the magnetization rotates 
along the easy axis inside the domain wall. One can cal- 
culate the change of tp across the domain wall to obtain 
Atp = \W h = (Vl + P' 2 - l)//3'- When /3' < 1, one 
finds that Atp — (3' /2 <C 1, which indicates that tp can 
be approximated as a constant within a domain wall as 
shown in Fig.[IJa), where the magetization always lies in 
the x-y plane [see also Fig.Q] (d)]. The solution, there- 
fore, returns to the Neel wall caseJ^ In contrast, when 
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ft' > 1, tp oscillates within a wall and the wall structure 
becomes more like a one-dimensional vortex with strong 
correlation between 9 and ip [see Fig.QJb) and (e)]. In fer- 
romagnetic metals, ft (~ 0.001-0.01) is smalliS^ and A° s 
is large due to the strong dipole-dipole interaction, hence 
the vertical spin stiffness is unimportant. However, in fer- 
romagnetic semiconductors, e.g., GaMnAs, A® s would be 
small 3 ii22 and ft can be large, ft' ~ ft ~ l) 20 i 23 therefore, 
the domain wall structure can behave like a vortex. 

When K± > 1 , the behavior of <p in a steady solution 
is determined by the competition between the hard-axis 
anisotropy, K±, and the transverse stiffness is propor- 
tional to ft'. Then ip tends to be locked to ip = nir (n 
is integer) and deviates from being a linear function of 
position, forming soliton-like structure or a staircase be- 
havior [Fig. [He) and (f)]. In a domain wall, modulation 
of 6 occurs when <j> is close to nir, but 8 tends to be close 
to n'n (n' is integer) when tp changes in order to lower the 
energy cost due to the hard axis anisotropy. Thus a do- 
main wall has an oscillating structure as seen in Fig.QJc) 
and (f ) . The number of oscillation increases as ft' is en- 
hanced. As far as we find numerically, the number of 
oscillation also depends on the initial condition of d x 9 
and d x tp at the boundary of the wall. Finally, we should 
point out that the domain- wall solution for K± ^ is the 
same as that for K± = in the absence of the vertical 
spin stiffness [see Fig.QJa) and (d)]32 

As is well-known^ the dynamics is strongly affected 
by the structure. For instance, vortex walls are easier 
to move than planar domain walls in the current-driven 
case, because of the perpendicular component near the 
vortex core. We may therefore expect that the wall 
for finite ft' > 1, e.g., in GaMnAs^ would have even 
lower threshold current due to the structure change aris- 
ing from the transverse exchange torque. Inclusion of 
the transverse stiffness in the micromagnetic simulations 
is thus crucially important in systems with strong spin- 
orbit interaction. 



V. SUMMARY 

In summary, we have derived the KSBEs in ferromag- 
netic systems based on the s-d model in the presence 
of the inhomogeneity of the magnetization. We analyti- 
cally solved the KSBEs and derived the spin torque due 
to the spin polarization of the itinerant electrons. The 
current-induced spin torque from the first-order magne- 
tization gradient is consistent with the previous works. 
We found that the second-order gradient of the magneti- 
zation inhomogeneity gives rise to an effective magnetic 
field that is perpendicular to the spin stiffness field. This 
vertical spin stiffness is proportional to the nonadiabatic 
parameter ft. We showed that the new term modifies the 
domain wall structure and causes magnetization rotation 
along the easy axis. The vertical spin stiffness is expected 
to be crucially important in ferromagnetic semiconduc- 
tors, and needs to be included in the LLG equation in 



numerical simulations on the magnetization dynamics. 
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Appendix A: Fourier transformation of Eq. ([7]) 

The Fourier transformations of the functions in Eq. ([7|) 
respect to the relative coordinate are calculated as fol- 
lows. For simplification, we omit all the temporal coor- 
dinates without leading to any ambiguity. 

The Green function at momentum k is defined as 

G(R,k) = J dre~ ikr G(R, r). (Al) 

Then, the electric potential energy term can be written 

as 

J dre- lk r [/ e (l)G(R, r) 

= j dre^ jk r e5 a Rf/ e (R)G(R,r) 

= e ihd k d«e J dke^ lk r [/ e (R)G(R, r) 

= e i i d ° a % e U e (R)G(R,k). (A2) 
Terms with Hamiltonian read 

J dre- 4k r i/ (Pi,ri)G(R,r) 

= j dre- lk - r e^ tfo(iP R + p,R)G(R,r) 
= e 4^9R C J dre- ik r iJ (iPR + P, R)G(R, r) 

= e <ia.8£° tf (ip R + k, R)G(R, k) 

- e ^(^ 9 R°-^ a "°)i7 (k,R)G(R,k), (A3) 

and 

j dre- jk r [G(R,r)ff (- P 2 ,r 2 )] 

= j dre- lk r e-^* [G(R,r)ffoH ^ + P> R )1 
= e -4^R° j dre- ik r [G(R,r) J ff (-i P R + P,R)] 

= e -^ 9kS - n G(R,k)tf (-i p R +k,R) 

= e -^( e ^R Q ^ 9 "°)G(R,k)7J (k,R). (A4) 



For the integral terms, one has Here, the time integral / 1H3 is omitted for simplification. 

Similarly, one can show 

J dre- lk r J dr 3 S(l, 3)G(3, 2) 



J dre- lkr J dr 3 e a * a ^E(Il,ri-r3) 



xe " 9 «G(R, r 3 - r 2 
= / dre 



/ dre"*- / dr*^* | ^^'^-^(R, k ') / d ^j *»E(1,3)G(3,2) 

xe ^a£ | *^ e ik''.( r3 -r 2 ) G(R) k // } = e^ 9 -~«>G(R, k )S(R, k). (A6) 



xe ik".(r 3 -r 3 )j G ( Rj k ') S ( R) k ") 



6 2 ( 9 k a £)£(R ; k)G(R, k). (A5) In these equations, we use the notation d R i9 k = V R • V k . 
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